Testing quantum nonlocality by generalized quasiprobability functions 
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We derive a Bell inequality based on a generalized quasiprobability function which is parameter- 
ized by one non-positive real value. Two types of known Bell inequalities formulated in terms of 
the Wigner- and Q-functions are included as limiting cases. We investigate violations of our Bell 
inequalities for single photon entangled states and two-mode squeezed vacuum states when vary- 
ing the detector efficiency. We show that the Bell inequality for the Q-function allows the lowest 
detection efficiency for violations of local realism. 
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I. INTRODUCTION 

Ever since the famous arguments of Einstein-Podolski- 
Rosen (EPR) [l[, quantum nonlocality has been a cen- 
tral issue for understanding the conceptual foundations 
of quantum mechanics. Quantum nonlocality can be 
demonstrated by the violation of Bell inequalities (Bis) 
which are obeyed by local-realistic (LR) theories. Re- 
alizations of BI tests are thus of great importance in test- 
ing the validity of quantum theories against LR theories. 
In addition, BI tests play a practical role in the detec- 
tion of entanglement which is one of the main resources 
for quantum information processing. Bell inequality tests 
for 2-dimensional systems have already been realized [|[ , 
while BI tests in higher-dimensional and continuous vari- 
able systems remain an active area of research [H, Q . 

Phase space representations are a convenient tool for 
the analysis of continuous variable states as they provide 
insights into the boundaries between quantum and clas- 
sical physics. Any quantum state p can be fully charac- 
terized by the quasiprobability function defined in phase 
space [f| . In contrast to the probability functions in clas- 
sical phase space the quasiprobability function is not al- 
ways positive. For example, the Wigner-function of the 
single photon state has negative values in certain regions 
of phase space Q • Since the negativity of the quasiprob- 
ability function inevitably reflects a non-classical fea- 
ture of quantum states, the relation between negativ- 
ity of quasi-probabilities and quantum nonlocality has 
been investigated [1, Q . Bell argued [1] that the original 
EPR state will not exhibit nonlocality since its Wigner- 
function is positive everywhere and hence serves as a clas- 
sical probability distribution for hidden variables. On the 
other hand, Banaszck and Wodkicwicz (BW) showed how 
to demonstrate quantum nonlocality using the Q- and 
Wigner-functions @. They suggested two distinct types 
of Bis, one of which is formulated via the Q-function and 



referred to in this paper as the BW-Q inequality while 
the other is formulated using the Wigner-function and 
is referred to as the BW-W inequality. Remarkably, the 
BW-W inequality was shown to be violated by the EPR 
state This indicates that there is no direct relation 
between the negativity of the Wigner-function and non- 
locality. 

Quasiprobability functions can be parameterized by 
one real parameter s 0, [Io| 



W(a;s) = 



7r(l 



-Tr[pU(a;s)}, 



(1) 



where U(a;s) = S^Lo(( s + l)/( s _ l)) n \ a i n )( a i n \, an d 
\a,n) is the number state displaced by the complex vari- 
able a in phase space. It is produced by applying the 
Glauber displacement operator D{a) to the number state 
\n). We call W (a; s) the s-parameterized quasiprobabil- 
ity function which becomes the P-function, the Wigner- 
function, and the Q-function when setting s = 1,0,-1 
[lfj, respectively. For non-positive s the function T4^(cv; s) 
can be written as a convolution of the Wigner-function 
and a Gaussian weight 



W(a; s) = 4, / d 2 (3 W(J3) exp 



(2) 
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This can be identified with a smoothed Wigner-function 
affected by noise which is modeled by Gaussian smooth- 
ing [ll|, [12, EH . Therefore decreasing s reduces the nega- 
tivity of the Wigner function and is thus often considered 
to be a loss of quantumness. For example, the Q-function 
(s = —1), which is positive everywhere in phase space, 
can be identified with the Wigner function smoothed over 
the area of measurement uncertainty. 

The purpose of this paper is to propose a method for 
testing quantum nonlocality using the s-parameterized 
quasiprobability function. We will firstly formulate a 
generalized BI in terms of the s-parameterized quasiprob- 
ability function in Sec. [Til This will lead us to a s- 
parameterized Bell inequality which includes the BW-Q 
and the BW-W inequalities as limiting cases. We will 
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then present a measurement scheme to test Bis using im- 
perfect detectors in Sec. IIHI The measured Bell expecta- 
tion value can be written as a function of the parameter 
s and the overall detector efficiency 77. In Sec. IIVI viola- 
tions of Bis will be demonstrated for single-photon entan- 
gled states and in Sec. [V] for two- mode squeezed vacuum 
states. We find the range of s and 77 which allows observ- 
ing non-local properties of these two types of states. We 
will show that the test involving the Q-function permits 
the lowest detector efficiency for observing violations of 
local realism. We also find that the degree of violation is 
irrespective of the negativity of the quasiprobabihty func- 
tion. Finally, in Sec. IVI1 we discuss the characteristics 
and applications of the s-parameterized BI. 



II. GENERALIZED BELL INEQUALITIES OF 
QUASIPROB ABILITY FUNCTIONS 

We begin by formulating a generalised BI in terms of 
quasiprobabihty functions. Suppose that two spatially 
separated parties, Alice and Bob, independently choose 
one of two observables, denoted by Ax, A 2 and Bx, B 2 
respectively. No restriction is placed on the number of 
possible measurement outcomes (which may be infinite). 
We assume that the measurement operators of the local 
observables Ax , A 2 , B\ , B 2 can be written as 

A a =d(a a ;s), B a = 6{(3 b ;s), for a, b =1,2 

using a Hermitian operator 



~ j (1 - s)fl(a; s) + st if-l<s<0, 

(a '' S) ~ I2n(a;s)-1 ifs<-l, 



(3) 



parameterized by a real non-positive number s and an 
arbitrary complex variable a. Here 1 is the identity op- 
erator. The possible measurement outcomes of 0(a; s) 
are given by its eigenvalues, 



(l-s)(f^) n + S 
2(f±lf-l 



if — 1 < s < 0, 
if s < -1, 



(4) 



and their eigenvectors are the displaced number states. 
The maximum and minimum measurement outcomes of 
0(a; s) for any non-positive s are A max = 1 and A m ; n = 
— 1, respectively. For s = we have 0(a; 0) = II(a; 0) = 
S^Lo( — l)"! ! n )( a J n L the displaced parity operator, 
while for s = — 1 we find that 0(a; — 1) = 2|a)(a| — 1 
projects onto the coherent states. 

A Bell operator can be constructed using the measure- 
ment operators A a , B b by way of a construction similar 
to the CHSH combination 



B = Cx,x + Ci, 



C 2 



c 2 



(5) 



where C a & = A a ® B^ is the correlation operator. 
Since the expectation values of the local observables are 



bounded by \{A a )\ < 1 an d | (■£?&} | < 1 for any non- 
positive s, the expectation value of the Bell operator 
defined in Eq. © is bounded by = \B\ < 2 in 

LR theories. Note that the expectation value of II(q;; s) 
for a given density operator p is proportional to the s- 
parameterized quasiprobabihty function 0, [Tpj 



W(a;s) 



-Tr[pfl(a; s)} 



tt(1 - s) 

v 1 n=0 v ' 

from which both the Wigner function and the Q-function 
can be recovered by setting s — and s = — 1, respec- 
tively. We do not consider the case s > when the 
eigenvalues of II(a; s) are not bounded. We thus obtain 
the following generalised BI 
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{-Ks<0} 



ir 2 (l-s) 4 



[W(ax,/3x;s)+W(ax,P2;s) 



+ W(a 2 ,/3x;s) - W(a 2 ,f3 2 ;s)} + tts(1 - s) 2 



x [W{ax;s)+W(/3x;s)]+2i 



<2, 



(7) 

\B\{.<-i} = \K 2 (l-s) 2 [W(ax,Px;s) + W(a 1 ,02;s) 
+ W(a 2 ,px;s)-W(a 2 ,l3 2 ;s)}-2Tr(l-s) 
x [W{ax;s) + W(f3x;s)] + 2\ < 2, 

where W(a, (3; s) = (4/tt 2 (1 - s) 2 ; )Tr [pTL(a; s) <x> n(/3; s)] 
is the two-mode s-parameterized quasiprobabihty func- 
tions, and W(a; s) and W(/3; s) are its marginal distribu- 
tions. We call Eq. (J7]) the s-parameterized Bell inequality 
for quasiprobabihty functions. This BI is equivalent to 
the BW-W inequality when s = which has the form of 
the standard CHSH inequality [3], and the BW-Q in- 
equality when s = — 1 in the form of the BI proposed 
by Clauser and Horn (CH) [l5|. In these cases the cor- 
responding generalized quasiprobabihty function reduces 
to the Wigner-function W(a,/3) = W(a,@;0) and the 
Q-function Q(a,/3) = W{a,f3;— 1), respectively @. 



III. TESTING QUANTUM NONLOCALITY 

In this section we present a scheme to test quantum 
nonlocality using the s-paramctcrized Bis. For a valid 
quantum nonlocality test the measured quantities should 
satisfy the LR conditions which are assumed when deriv- 
ing Bis. Thus we here employ the direct measurement 
scheme of quasiprobabihty functions using photon num- 
ber detectors proposed in [l2l ], 

A pair of entangled states generated from a source of 
correlated photons is distributed between Alice and Bob, 
each of whom make a local measurement by way of an 
unbalanced homodyne detection (see Fig. [1]). Each lo- 
cal measurement is carried out using a photon number 
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FIG. 1: The optical setup for the BI test. Each local mea- 
surement is carried out after mixing the incoming field with 
a coherent state (denoted by |£) for Alice and \S) for Bob) in 
a beam splitter (BS) of high transmissivity T. The photon 
number detectors (PNDs) have efficiency 7)4. 



ence or absence) measurements (s = —1). 

Let us now consider the effects of the detector efficien- 
cies T). If the true photon number distribution is given by 
P(n), then the measured distribution can be written as 
a function of the overall detection efficiency 77 = rj^T as 
PvH = En= m Pn(^)(l r V y-™ V m M. For a = 
the measured quasiprobability function is 



W„(0;a) 



8 + 1 



7r(l — s) ^— ' V s — 1 

v ' m=0 v 



E 



7r(l - s) 



n=0 



J2 ( 1 -v+v 



P v (m) 

ill 

s-1 



P(n) 



W 



_ W(0;s') 
V 



(8) 



detector with quantum efficiency rjd preceded by a beam 
splitter with transmissivity T. Coherent fields |£) and \S) 
enter through the other input ports of each beam split- 
ter. For high transmissivity T — > 1 and strong coherent 
fields £,<5 — ► 00, the beam splitters of Alice and Bob 
can be described by the displacement operators D{a) 
and D{(3) respectively, where a 
/3 = Sy(l-T)/T^ 
with n — ^2 n n\n)(n\ the photon number operator are 
performed on the outgoing modes using perfect photon 
number detectors. Then the expectation value directly 
yields the value of the s-paramctcrized quasiprobability 
function at the point in phase space specified by the com- 
plex variables a and /3. For example, the Wigner func- 
tion can be obtained by the parity measurements (— 1)™ 
(s = 0) and the Q-function by on-off (i.e. photon pres- 



e^/CT- T)/T and 
Measurements ((s + l)/(s — l)) r 



The s-parameterized quasiprobability function measured 
by a detector with efficiency 77 can therefore be identi- 
fied with the quasiprobability function with parameter 
s' = — (1 — s — r])/r]. Other sources of noise (e.g. dark 
counts and mode mismatch) could be included into this 
approach but are neglected here for simplicity. 

Finally, the expectation value of observable Q is given 

as 



(6(a;s)) v = 



7T(1-S) 



-W(a;s') + s if-l<s<0, 
W(a;s')-l if s < -1, 



(9) 



where (•),, represents the expectation value obtained by 
measurement with efficiency rj. Note that (0) is the sta- 
tistical average of directly measured data without posts- 
election. The expectation value of the Bell operator ([5]) 
written as a function of s and 77 is given by 



J 



(S{-Ki<0})ti — 



(&{s<-i})v - 



n 2 (\-s) A 



W(ai,/3i 



I — S — T] 

V 

V 



TTS(I-S) 2 
+- " - 

V 



1 — s — rj 



)+W{ct 2 ,fc- 



1 — s — rj . 



Wiav,- 1 S V ) + W((3v,- 



V 

1 — S — 7] , 



T) 



W{ai,{h;—- — - — -) + W(ai,{3 2 ;-- — - — -) + W(a2,Pv,- 



V 

1 - s - T] , 
V 



1 - s - 77 
V 



2tt(1 - s) 

T) 



+ 2s 2 , 



(10) 



Note that the Bell expectation values in Eq. (|10j) for s = 
and s = — 1 give the same results as tests of the BW-W 
and BW-Q inequalities, respectively. 



IV. VIOLATION BY SINGLE PHOTON 
ENTANGLED STATES 

We investigate violations of the s-parameterized BI (0 
for the single photon entangled state [I?) 



I*) 



V2 



(|0,1> + |1,0)), 



(11) 
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FIG. 2: Maximum Bell expectation value \B\ = \{B)\ for the 
single photon entangled state. Only the range of parameters 
s and detector efficiencies r\ with \B\ > 2 is shown. 



where \n, m) is the state with n photons in Alice's mode 
and m photons in Bob's mode. This state is created by a 
single photon incident on a 50:50 beam splitter. Its two- 
mode s-parametcrized quasiprobability function is given 
by 



Wy(a,/3;s) 



1 



7T 2 (1- S ) 2 

x exp 



1-s (1-s) 2 
2(M 2 + I/3| 2 ) 



(12) 



and its marginal single-mode distribution is 

Wv(a;s) = (l/7r)(2-2?7 + V|a| 2 )exp[-2?7|a| 2 ]. (13) 

Note that for > s > — 1 Eq. (fT2")l has negative values in 
certain regions of phase space but for s = — 1 it becomes 
the Q-function Wq,(a,(3;— 1) > 0. 

The maximum expectation values |£>| max = |(£>)| max 
are obtained for properly chosen on, ai, (3\, fii- Figure [2] 
shows the range of parameters s and detector efficiencies 
rj for which the BI is violated, |i3| max > 2. Interestingly 
the degree of violation is not directly related to the neg- 
ativity of the quasiprobability functions. The test of the 
BI using the Q-function (s = —1) yields strong viola- 
tions and is most robust to detector inefficiencies. This 
is because the observable ([3]) becomes dichotomized at 
s = — 1 corresponding to detection of none vs. some pho- 
tons. For a given s, the amount of violation decreases 
with decreasing rj. The minimum value of r\ indicates 
the required detector efficiency for a successful nonlocal- 
ity test [l8j |. For example, the minimum bound is about 
83% for the Q-function (s = —1). We also find the mini- 
mum parameter s which allows demonstrating quantum 
nonlocality for a given detector efficiency. For example, 
for a perfect detector (rj = 1), the corresponding BI is 
violated when s > —1.43. 



VIOLATION BY TWO-MODE SQUEEZED 
STATES 



We consider the two-mode squeezed vacuum states 
(TMSSs), i.e. a continuous variable entangled state writ- 
ten as 



|TMSS) 



sech r tanh" r\n, n), 

n=0 



(14) 



where r > is the squeezing parameter. It can be re- 
alized for instance by non-degenerate optical parametric 
amplifiers fl9j | . In the infinite squeezing limit r — > oo, the 
TMSS becomes the normalized EPR state which is the 
maximally entangled state associated with position and 
momentum 

For a non-positive s the quasiprobability function of 
the TMSS is given by 



W T MSs(a,f3;s) 



exp 



tt 2 R(s) V R ( s 
- smh2r(a(3 + a* P*)} 



{S(s)(\c 



and its marginal single-mode distribution is 



IT- 



TMSS 



a;s) 



exp 



S(s) 



-\P\ 2 ) 
(15) 

(16) 



where R(s) = s 2 — 2s cosh 2r + 1 and S(s) — cosh 2r — s. 
Note that these are positive everywhere in phase space. 
In Fig. 03a) violations of the s-parameterized BI are 
shown for TMSSs. The test using the Q-function (s = 
— 1) is most robust with respect to detector inefficien- 
cies. The amount of violation shows different tendencies 
depending on the squeezing parameter r. In the case of 
low squeezing rates, i.e. when the amplitudes of small-n 
number states are dominant, the violation is maximal if 
we choose the Q-function (s = —1) as shown in Fig.[3](b). 
This implies that the dominant contribution to the vio- 
lation comes from correlations between the vacuum and 
photons being present. For larger squeezing rates r > 1.2, 
the violation reaches a maximal value B ~ 2.32 when we 
test the Wigner-function (s = 0) [2(| ■ This indicates 
that the parity measurements are effective for verifying 
higher-order number correlations. However, the parity 
measurements require very high detector efficiency as 
shown in Fig. [3jb). The range of s within which one 
can demonstrate nonlocality becomes narrower around 
s = and s = — 1 with increasing squeezing rate r. This 
is because the observable ^ is dichotomized at s = 
and s = — 1. 



VI. DISCUSSION AND CONCLUSIONS 

We demonstrated that quantum nonlocality has no 
direct relation to the negativity of ,s-parameterized 
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FIG. 3: Demonstration of quantum nonlocality for TMSSs. (a) Maximum Bell values are shown for different squeezing r in 
the range of s and rj where the BI is violated, (b) Violation of the BI as a function of the squeezing r for different s and r\ = 1 
(solid line), rj — 0.95 (dashed line), and rj — 0.9 (dotted line). 



quasiprobability functions. In fact the Q-function (s — 
— 1) which never becomes negative can still be used to 
verify non-local properties as we showed in Fig. [2] and 
yields strong violations of the corresponding BI. This im- 
plies that the quantum properties of nonlocality and neg- 
ativity of the quasiprobability functions should be con- 
sidered distinct features of quantum mechanics. Further- 
more we showed that the Q-function test allows the low- 
est detector efficiency for demonstrating quantum non- 
locality. For example, it requires only rj sa 83% for a 
single photon entangled state and rj « 75% for TMSSs 
with r = 0.4 to detect nonlocality. This indicates that 
two-mode correlations between vacuum and many pho- 
tons can be more robust to detector inefficiencies than 
correlations between vacuum and a single photon. 

The parameter s determines the characteristics of the 
detected non-local correlations. For example, if we 
choose s = — 1 the violation of the BI exhibits only cor- 
relations between vacuum and photons. In order to test 
higher-order photon number correlations we need to in- 
crease s to zero, so that the factor ((s + l)/(s — 1))" 
multiplied to the photon number probability increases in 
Eq. ©. Although parity measurements (s = 0) allow to 
detect higher-order correlations effectively, they also re- 
quire very high detector efficiencies as shown in Fig. [3l If 
we properly choose a certain parameter — 1 < s < 0, e.g. 
s = —0.7, we can detect higher-order correlations with 
a lower detector efficiency than that required for test- 



ing the BI using the Wigner function. However, we note 
that the violation of the BI with s = —0.7 disappears 
with increasing squeezing rate as shown in Fig. [3f b) ; this 
restricts the possible applications to schemes using light 
that contains only a few photons. 

Let us finally discuss whether we can regard de- 
coherence effects as changes to s. Interactions with 
the environment and detection noise tend to smoothen 
quasiprobability functions. For example, when solv- 
ing the Fokker-Planck equation for the evolution of the 
Wigner-function of a system interacting with a thermal 
environment one obtains plj 

W(a,r) = I ^J dW th (/3)W^^^,T = o)(17) 

Here the parameters r(r) = y/1 — e fT and t(r) = 
VfT^ are given in terms of the energy decay rate 7, 
and 

7r(l + 2n) \ 1 + 2n J 

is the Wigner function for the thermal state of average 
thermal photon number n. The effect of the thermal 
environment is then identified with temporal changes of 
the parameter 
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Therefore one might be tempted to consider an envi- 
ronment in a thermal state as giving rise to a temporal 
change in s in Eq. (|3|). However this idea is not applica- 
ble to tests of quantum nonlocality. The s-parameterized 
BI is derived for observables ([3]) which contain s as a de- 
terministic value of LR theories. Thus the local- realistic 
bound is no longer valid when dynamical observables are 
considered (even though they give the same statistical av- 
erage). However, this idea might be useful for witnessing 
entanglement [23 ], 

In summary, we have formulated a BI in terms of the 
generalized quasiprobability function. This BI is param- 
eterized by a non-positive value s and includes previously 
proposed Bis such as the BW-W (s = 0) and the BW-Q 
(s = —1) inequalities Q. We employed a direct measure- 



ment scheme for quasiprobability functions [12[ to test 
quantum nonlocality. The violation of Bis was demon- 
strated for two types of entangled states, single photon 
entangled and two-mode squeezed vacuum states. We 
found the range of s and 77 which allow the observation 
of quantum non-local properties. We discussed the types 



of correlations and their robustness to detection ineffi- 
ciencies for different values of s. We also demonstrated 
that the negativity of the quasiprobability function is 
not directly related to the violation of Bis. The realiza- 
tion of s-parameterized BI tests is expected along with 
the progress of photon detection technologies [23| in the 
near future. Our investigations can readily be extended 
to other types of states like photon subtracted gaussian 
states [24], [25| , or optical Schrodinger cat states |2fi| . 
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